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We present a theoretical framework which is generally ap- 
plicable to the study of time scales of activated processes in 
systems with Brownian type dynamics. This framework is ap- 
plied to a prototype system: magnetization reversal times in 
the 2D Ising model. Direct simulation results for the magneti- 
zation reversal times, spanning more than five orders of mag- 
nitude, are compared with theoretical predictions; the two 
agree in most cases within 20%. 
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I. INTRODUCTION 

Activated processes that can be described with some 
type of Brownian dynamics are abundant in the world 
around us. Well-known examples are the nucleation of 
droplets in an undercooled gas or of crystals in an un- 
dercooled liquid, chemical reactions and the escape of 
a protein from a misfolded state. A prototype system 
to study such phenomena numerically is the well-known 
Ising model. Above the so-called critical temperature, 
in absence of an external magnetic field, up- and down- 
pointing spins are roughly equally abundant. Below the 
critical temperature, the system prefers to be in either 
of two states: one state with a positive magnetization in 
which most spins are pointing up, and one state with a 
negative magnetization. As long as the system size re- 
mains finite, reversals of the magnetization - transitions 
between positive and negative magnetization - are possi- 
ble and will occur at a certain average frequency. These 
processes are activated, since configurations with mag- 
netization close to zero have a higher free energy than 
typical configurations with a magnetization close to ei- 
ther of the equilibrium values. 

In this manuscript, we study the time scales associ- 
ated with magnetization reversal. A theoretical frame- 
work is outlined which is generally applicable to activated 
processes in systems with Brownian type dynamics, and 
compared to high-accuracy computer simulations. From 
a practical point of view magnetization reversals are also 
of great interest because of applications in memory de- 
vices and the like. One wants to have rapid switching 
of magnetization under reversals of an external field, but 
no spontaneus reversals of the magnetization, even if the 
external field has been turned off. In the literature much 
attention has been paid to reversal time distributions in 
the presence of a driving field [1,2], but spontaneous re- 



versals in the absence of a field have hardly been studied. 
Here we consider the latter case for the prototypical case 
of an Ising model with periodic boundary conditions. We 
identify the leading scenario for reversals of the magneti- 
zation and show that the process may be described to a 
good approximation by a one-dimensional diffusion pro- 
cess over a potential barrier. 

Our manuscript is organized as follows. In section II, 
we describe the model that we study in detail. Next in 
section III, we outline the theoretical framework which 
is generally applicable to activated processes in systems 
with Brownian type dynamics. We then apply this frame- 
work to our prototypical model — magnetization reversal 
in the Ising model. In section IV we compare the theo- 
retical predictions with high-accuracy computational re- 
sults. 



II. DETAILED DESCRIPTION OF THE MODEL 

We consider the Ising model on a B x L rectangular 
lattice with periodic (helical) boundary conditions, with 
the Hamiltonian 



H 



(id) 



(1) 



in which Oj = ±1 is the spin at site i and J is the coupling 
constant. The summation runs over all pairs of nearest- 
neighbor sites; those of site i are j = i ± 1 modulo N 
and j = i ± B modulo N, with N = BL. The mag- 
netization is defined as M = X^^ii ^ can take values 
M = —N, —N + 2, . . . , N; all through this manuscript, 
we restrict ourselves to systems in which both B and L 
are even. As a consequence, M takes only even values, 
and summations over a range of possible magnetizations 
only run over even numbers, with an increment of 2. 

A well-known property of the model is that at suffi- 
ciently low temperatures, the distribution of the magne- 
tization M = '^2 i (Ji becomes bimodal: the spins in a con- 
figuration tend to align around the two preferred values 
±Mo. For infinite systems this occurs at all temperatures 
below the critical temperature T c = J/(0.44069/cb), for 
finite systems this range starts at slightly different tem- 
peratures. 

The system evolves in time according to single- spin- flip 
dynamics with Metropolis acceptance probabilities [4] . If 
Cj is the configuration after i proposed spin flips, a trial 
configuration C' i+1 is generated by flipping a single spin 
at a random site. This trial configuration is then cither 



1 



accepted (Cj+i = C- +1 ) or rejected (Cj+i = Cj); the 



acceptance probability is given by 

P a = min [l,exp(-/3(£;(^ +1 ) - E(C\))] , 



(2) 



in which (3 = l/(fcsT) with Boltzmann constant fee and 
temperature T. The time scale is set such that in one unit 
of time, on average each spin is proposed to be flipped 
once. So in our system, in one unit of time we perform 
BL Monte Carlo steps. 




FIG. 1. Snapshots of the the transition from a state with 
most spins down to a state with most spins up in the 16 x 64 
Ising model, taken at equal time intervals of 500 attempted 
spin flips per site, at inverse temperature f3J = 0.5. 

Our current interest in the Ising model with single- 
spin-fhp dynamics stems from the fact that in finite sys- 
tems, the configurations will occasionally switch between 
states in which the magnetization is either negative or 
positive, through an activated process. The dominant 



pathway at low temperatures consists of the formation 
of a single pair of closed interfaces in the shorter peri- 
odic direction (for B ^ L), which perform a relative dif- 
fusive motion around the longer periodic direction and 
annihilate after meeting each other through the periodic 
boundary. A series of snapshots, illustrating this process, 
is presented in figure 1. 



III. THEORETICAL FRAMEWORK 

To study the behavior of times between magnetization 
reversals at temperatures below the critical one we may 
consider an ensemble of a large number of systems pre- 
pared in states with a magnetization close to the equi- 
librium value —Mq and study the rate at which these 
systems reach the value M = M , after which they are 
removed from the ensemble. Alternatively, the behavior 
of times between zero-crossings of the magnetization can 
be studied, as discussed below. The spin-flip dynamics 
described above may be represented by a master equa- 
tion for the probability distribution -P(cr) of finding a 
system in the state a at time t. Due to the huge number 
of possible states this master equation cannot be solved 
analytically or even numerically for system sizes of prac- 
tical interest. Therefore, as an approximation we assume 
that we may replace the exact master equation by an ap- 
proximate master equation for the probability P(M, t) of 
finding a system with magnetization M at time t. The 
form of this equation is 



dP(M,t) 
dt 



— Fm,M+2P(M + 2, t) + Tm,M-2 P(M-2,t) 



— (Tm+2,M + ^M-2,M)P(M,t), 



(3) 



with Tm',m the transition rate from M to M' . It ignores 
the fact that the actual transition rates will depend on 
the geometry of the state under consideration, with the 
idea that in typical cases spin states of very similar ge- 
ometry will dominate the set of states with given M. 
Besides on M and (3 the transition rates will depend on 
the geometry of the system, that is on B and L. To esti- 
mate the scaling behavior of these dependencies we make 
two further simplifying assumptions: 

• First of all, requiring B < L, we assume that states 
with M ±Mo consist of a single strip of opposite 
magnetization separated by two phase boundaries 
of length B from the majority spin phase. The 
relevant changes of M then will be caused by dis- 
placements of these boundaries due to flips of spins 
along them. The total number of spins available 
for this will be proportional to B, hence the tran- 
sition rates should also be proportional to B and 
independent of L. 

This approximation will not be valid for values of 
M that are too close to ±M , as for these the op- 
posite magnetization will typically be found in a 
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closed cluster rather than in a strip. It will turn out 
though that the contributions from these M- values 
to the reversal frequency are very small for systems 
of reasonable length. Therefore using the approxi- 
mation r ~ B also here does not harm. Further we 
neglect the possibility of having more than one strip 
of opposite magnetization. When the temperature 
gets close to the critical one, when B becomes small 
or when L becomes very large, this may not be a 
good approximation. We will come back to this in 
section V. In our further theoretical treatment we 
will assume that we are in a situation where the 
single-strip approximation is justified. 

• Further we neglect changes due to fluctuations 
of the magnetization caused by the growth and 
shrinkage of small clusters of the minority spin 
type. On average these do not contribute to magne- 
tization reversals, so this approximation should be 
allowed. In case one considers the first passage fre- 
quency through zero magnetization there is a small 
effect, to which we will come back in section V too. 

In order that the equilibrium distribution be a sta- 
tionary solution of the master equation we impose the 
condition of detailed balance 



M.M+2 
^M+2,M 



= cxp [P(F(M) - F(M + 2))] , 



(4) 



where (3F(M) = -InP e ,(M), with P eq (M) the equilib- 
rium probability of magnetization M. In section IV B it 
will be explained how this condition may be combined 
with simulation results for interface diffusion to obtain 
approximations for all transition rates that lead to a good 
overall prediction of the reversal times. 

The long-time reversal frequency as predicted by the 
master equation (3) follows as the largest eigenvalue of 
this equation, supplemented with an absorbing boundary 
at M = A, with either A = Mq or A = 0. The first choice 
corresponds to a real reversal of the magnetization, the 
second one to a first return to the value M = 0. After 
this the system will have equal probabilities to actually 
reverse its magnetization or to return to the equilibrium 
magnetization value it came from, hence this return fre- 
quency should be twice the reversal frequency. The ab- 
sorbing boundary condition is implemented by setting 
T*A—2,A equal to zero. 

The largest eigenvalue — v of Tm,m> m ec l- as weu 
as the corresponding eigenvector P (M), may be found 
by requiring that the net current away from magnetiza- 
tion M assumes the value vPq(M). Using conservation of 
probability one easily checks that this may be expressed 
as 

r M +2,M-Po(M) - T m ,m+2Po(M + 2) = v J2 P oM, (5) 
T A , A _ 2 P {A - 2) 



i<M 



In eq. (5) Po(m) on the right hand side may be approx- 
imated, up to a normalization factor, by exp(— (3F(m)), 
because the sum is dominated by the terms with small 
m-values, for which this approximation is excellent. This 
one may check in hindsight against the solution ob- 
tained. With this approximation the equation may be 
solved recursively for Pq{M) in terms of Po(A — 2) for 
M = A - 4, A - 6, with the result 



P (M) 



E 



A-2 



M<m<A-2 



- m+2,m 



exp [(3{F{m) — F{M))] 



E„ £ „e,p[-m„)i 
E„'<A-2 ex P [-WW)] 

The summations over m and M are dominated by values 
of m close to zero, combined with values for M close to 
—Mq, for which the sum over n is basically independent 
of m. Therefore the reversal frequency may be obtained 
as 

Ta,a-2Pq(A-2) 
Em<a-2Po(M) 

A-2 



\ - cxp[/3F(m)] x - 
2^ fi 2^ exp[-#F(n)] 



. r m _|_2.m 

. m— — Mq n<A — 2 



• (8) 



jm<A-2 



P (m) 



(6) 



The restriction of the summation over m to values 
m > —Mq is needed to avoid large spurious contributions 
from m < —Mq. The result in eq. (8) is well-known. It is 
usually derived by considering a state with a stationary 
current in which mass is inserted at a constant rate on 
one side (e.g. at M = —B x L in our case) and taken out 
as soon as it reaches the absorbing boundary (see e.g. [6], 
section IV E). In that case the replacement of the sum 
over n by a constant is exact. 



IV. SIMULATIONS AND RESULTS 

In order to apply the above theoretical framework to 
magnetization reversal times in the Ising model, the two 
ingredients required are: (i) the equilibrium probability 
P eq (M) to find the system in a state with magnetization 
M; and (ii) the transition rates Tm'.m from magnetiza- 
tion M to M 1 . We obtain these two ingredients via two 
different computational approaches. 



A. free energy landscape 

For various values of L and B and various tempera- 
tures we make histograms of the distribution of the mag- 
netization M = <7j and the energy. Since the prob- 
ability that a state with energy E occurs at an inverse 
temperature (3 is proportional to the Boltzmann weight 
exp(— (3E), a histogram made at a certain temperature 



3 



provides information about the probability distribution 
at nearby temperatures as well. We use the multiple his- 
togram method [5] to combine the information from sim- 
ulations at various temperatures. In this way, histograms 
for the magnetization over a wide range of temperatures 
can be obtained. As defined above, the free energy F(M) 
is related to the probability that a certain magnetiza- 
tion M occurs by P eq (M) = e -P F i M ) . The method of 
determining free energies of configurations constrained 
to some value of a coordinate along a pathway through 
phase space has been used extensively, for instance by 
Auer and Frcnkcl [7] . In figure 2 the free energy as a func- 
tion of the magnetization per spin m = M/N is plotted 
at two different temperatures and three different values 
of L/B. 
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FIG. 2. Free energy as a function of magnetization at in- 
verse temperatures /3J = 0.45 and /3J = 0.50. The number of 
spins is kept constant at 1024. The solid lines represent the 
32 x 32 system, the dotted lines represent the 16 x 64 system, 
and the dash-dotted lines represent the 8 x 128 system. The 
free energies at M = corresponding to equation (10) are 
also indicated. A constant is added to the free energy curves, 
such that the minimal value of F(M) is zero in all cases. 

In order to compare these results with theory we note 
that the probability P eq (0) of finding the system in a state 
with M — can be obtained from eq. (4.23) of chapter 
V of [8] as follows. First note that this equation gives 
the following approximation for the partition function 
of a periodic anti-ferromagnetic system with one phase 
boundary [9] , which is equivalent to a ferromagnetic sys- 
tem with an anti-periodic boundary: 



Zi 




A(k) - ^/A 2 (k)-c 2 



B/2 



-{ \A(k) + ^A 2 {k)-c 2 ^ 



(9) 



where 



A(k) = (1 + z 2 ) 2 - 2z(l - z' z ) cos 



irk 



c = 2z(l - z 2 ), 
z = tanh/3J, 

and the partition function for a homogeneous system 
(cf. (4.20) of [8]). Neglecting the possible interaction be- 
tween two phase boundaries (or interfaces) we conclude 
that the partition function Z 2 for a system with two in- 
terfaces is given by Z 2 /Z = (l/2)(Z 1 /Z a ) 2 . The factor 
1/2 arises because otherwise we would count each con- 
figuration twice: interchanging the locations of the two 
interfaces does not give a different configuration. If we 
leave the positions of the two interfaces completely free 
the magnetization is distributed almost uniformly over all 
possible even values between — M and M . Therefore, 
fixing M to leads to a reduction of the partition func- 
tion by a factor of M . Thus we arrive at the following 
result for P eq (0): 



M Z 
1 

8Mo 



( A(k) - ^A 2 (k)-c 2 

fe= 



B/2 



^ \A(k) + y/A*(k)-<? 



(10) 



For square systems this result has to be modified, be- 
cause there are two ways to make a strip with opposite 
magnetization: the interfaces can lie in the horizontal as 
well as in the vertical direction. This gives an additional 
factor 2 in the equation for P eg (0). 



B. interface diffusion coefficient 

The second ingredient for the theoretical framework in 
section III consists of the transition rates Tm'.m from 
magnetization M to M' . As discussed there, we are ac- 
tually mostly interested in the contribution to the tran- 
sition rates that arises from the diffusion parallel to 
the longer periodic direction of interfaces that span the 
shorter one. To estimate this diffusion coefficient nu- 
merically, we study systems with anti-periodic boundary 
conditions: the spins in two neighboring sites i < N — 1 
and j = i + 1 , or in i < N — B and j = i + B arc aligned 
if Ui = <jj, whereas the spins in two neighboring sites 
i = N - 1 and j = 0, or in i > N - B and j =i + B-N, 
are aligned if <r, = — aj. 

The interface location x(0) is initially defined as the 
magnetization M(0). As long as M £ [-0.8N,0.8N], 
steps Ax in the interface location are equal to changes 
AM in the magnetization, i.e., x(t+ At) —x(t) = M(t + 
At) — M(t). Once the interface gets close to the anti- 
periodic boundary, the magnetization does not uniquely 
determine the interface location. As soon as M is no 
longer in the interval [— 0.8iV, 0.8AT], we shift the anti- 
periodic boundary away from the interface over half the 
system size (which will bring the magnetization back in 
this range) and then continue. 
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In practice, we achieve this by switching from moni- 



toring M to monitoring M 1 = J2i=o~ a i ~ J2i=N/2 u ^ 
and \x(t + At) - x(t)\ = \M'(t + At) - M'(t)\. The sign 
of the steps in x(t) depends on the location of the in- 
terface: in the lower half this sign will be unchanged, 
whereas in the upper half it will be reversed. As soon as 
M' leaves the interval [-0.8N, 0.8N] we switch back to 
measuring steps in the original magnetization M. Also 
here, the sign of the steps in x(t) depends on the location 
of the interface. Note that x is not confined to the inter- 
val [—A, N]. The diffusion coefficient D is then obtained 
from the time-dependent interface location x(t) as 



D 



lim 

t— >oo 



((x(t)-x(0)f 



2t 



(11) 



We expect that the interface diffusion is independent of 
L, and that it grows linearly with B since the number of 
sites where a spin flip moves the interface grows linearly 
with its width; however, for small B corrections arise due 
to the periodicity (or helicity) of the boundaries. Conse- 
quently, we expect 



D(B,L,/3J) =g((3J)B + c. 



(12) 



For temperatures close to T c and small B, the system 
might occasionally contain more than a single interface. 
In that case, equation (11) would overestimate the diffu- 
sion coefficient of an interface. Since the additional free 
energy cost of interfaces increases linearly with B, this 
unwanted contribution to D decreases exponentially with 
B. Taking this into account, our approach is to measure 
for various values of (3J, B, and L the diffusion coeffi- 
cient D via eq. (11). Next we determine the function 
g(PJ) in equation (12); the results are plotted in figure 
3. Indeed we find that for increasing B and at tempera- 
tures not too close to T c the diffusion coefficient rapidly 
becomes independent of L within the ranges of L-values 
we considered. In our theoretical framework we then use 
as approximation for the jump rates 



min[rM,Af+2, Tm+2,m] — 



g(0J)B 



(13) 



This may be understood from the observations that on 
the one hand one has D = 4T, because the jumps in 
magnetization go by units of 2, on the other hand D 
satisfies (12) with B replaced by 2B because there are two 
interfaces. The above equation, in combination with eq. 
(4) and the free energy as a function of magnetization, 
specifies all the transition rates. Using (8) we can now 
predict the magnetization reversal times. The results are 
shown in Table I. 
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FIG. 3. Monte carlo measurements of the diffusion coef- 
ficient per interface length g(/3J), as a function of inverse 
temperature f3J . 



TABLE I. Estimated values for the magnetization reversal 



times. 


PJ 


32 x 32 


16 x 64 


16 x 32 


0.46 


1.66 


10 b 


1.74 


10 4 


1.19 • 10 4 


0.47 


1.12 


10 6 


4.21 


10 4 


3.15 ■ 10 4 


0.48 


8.88 


10 6 


1.25 


10 5 


9.96 • 10 4 


0.49 


7.38 


10 7 


3.97 


10 5 


3.31 • 10 5 


0.50 


6.49 


10 8 


1.31 


10 6 


1.12 • 10 6 


0.51 


5.83 


10 9 


4.19 


10 6 


3.61 • 10 6 


0.52 


5.31 


10 10 


1.23 


10 7 


1.05 • 10 7 


0.53 


5.77 


10 11 


3.78 


10 7 


3.11 • 10 7 


0.54 


7.69 


10 12 


1.23 


10 8 


9.60 • 10 7 


pj 


16 x 16 


8 x 64 


8 x 32 


0.46 


4.43 


10 3 


2.76 


10 3 


1.24 • 10 3 


0.47 


1.06 


10 4 


3.26 


10 3 


1.73 ■ 10 3 


0.48 


3.09 


10 4 


4.72 


10 3 


2.88 • 10 3 


0.49 


9.56 


10 4 


7.45 


10 3 


5.11 • 10 3 


0.50 


3.09 


10 5 


1.26 


10 4 


9.47 • 10 3 


0.51 


9.60 


10 5 


2.12 


10 4 


1.69 • 10 4 


0.52 


2.73 


10 6 


3.32 


10 4 


2.77 • 10 4 


0.53 


8.06 


10 6 


5.50 


10 4 


4.71 • 10 4 


0.54 


2.50 


10 7 


9.69 


10 4 


8.43 ■ 10 4 


0.55 


7.64 


10 7 


1.70 


10 5 


1.49 • 10 5 


0.56 


2.30 


10 8 


2.96 


10 5 


2.61 • 10 5 


0.57 


7.20 


10 8 


5.38 


10 5 


4.75 • 10 5 


0.58 


1.92 


10 9 


8.37 


10 5 


7.39 • 10 5 


0.59 


5.54 


10 9 


1.42 


10 6 


1.25 • 10 6 


0.60 


1.77 


io 10 


2.66 


10 6 


2.34 • 10 6 



•5 



C. magnetization reversal times 

We measure the magnetization every TV attempted spin 
flips. We look for events where the magnetization crosses 
from positive (more than half the spins up) to negative 
or vice versa between two consecutive measurements. 

We make histograms of the times between two occur- 
rences of magnetization reversal events. The histogram 
obtained in a system containing 16 x 64 spins, at the 
temperature corresponding to [3J = 0.45, is shown in 
figure 4. For comparison we also show a histogram of 
the times measured between the first time the system 
reaches a free energy minimum, and the first time after 
this it reaches the other minimum. We will come back to 
this in our discussion. 
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FIG. 4. Histogram of the times between zero-crossings of 
the magnetization (squares), and the times between first oc- 
currence of a free energy minimum (circles) in the 16 x 64 
Ising model at (3 J — 0.45. The solid lines depict the fit to the 
data, obtained as described below. 

The figure shows that at long times the decay function 
f(t) behaves as f(t) ~ exp(— t/r). Here, we focus on 
the long-time behavior, and are specifically interested in 
the escape time r. We obtain this quantity via a fitting 
procedure, in which we ignore the data up to a time to, 
chosen such that f(t) shows exponential time behavior 
for t > to- Then we determine the time t' at which half 
of the remaining events have taken place. The escape 
time t is then obtained from r = (t' — t )/ m(2). Instead 
we could have made a linear fit of all the data points 
beyond to, but this makes no significant difference. 

The resulting reversal times, for several system sizes 
and inverse temperatures, are presented in Table II and 
plotted in figure 5. 




0.5 0.55 
pj 



FIG. 5. Magnetization reversal times r as a function of in- 
verse temperature (3 J, for various system sizes. The diamonds 
represent systems with L = 16, the circles represent systems 
with L = 32, and the squares represent systems with L = 64. 



TABLE II. Directly measured values for the magnetization 
reversal times. 



pj 


32 x 32 




16 x 64 


16 x 32 




0.46 


2.09(6) 


• lO 5 


1.52(3) 


• 10 4 


1.37(1) 


• 10 4 


0.47 


2.02(6) 


• 10 6 


5.14(5) 


• 10 4 


4.18(4) 


• 10 4 


0.48 


1.97(6) 


•10 7 


1.68(1) 


• 10 5 


1.32(1) 


■ 10 5 


0.49 






5.20(5) 


• 10 5 


4.26(6) 


• 10 5 


0.50 






1.62(5) 


• 10 6 


1.33(2) 


• 10 6 


0.51 






5.2(1) ■ 


10 6 


4.28(6) 


• 10 6 


0.52 






1.59(5) 


•10 7 


1.32(3) 


•10 7 


0.53 






4.5(2) ■ 


10 7 


4.0(1) • 


10 7 


0.54 






1.40(6) 


• 10 8 


1.13(4) 


• 10 8 




16 x 16 




8 x 64 




8 x 32 




0.46 


6.83(7) 


• 10 3 


1.21(1) 


• 10 3 


1.23(1) 


• 10 3 


0.47 


1.81(2) 


■ 10 4 


2.10(2) 


• 10 3 


2.08(2) 


• 10 3 


0.48 


4.8(2) • 


10 4 


3.75(3) 


■ 10 3 


3.61(4) 


• 10 3 


0.49 


1.35(4) 


• 10 5 


6.70(7) 


■ 10 3 


6.30(6) 


• 10 3 


0.50 


4.45(4) 


• 10 5 


1.20(1) 


• 10 4 


1.11(1) 


• 10 4 


0.51 


1.26(1) 


• 10 6 


2.07(2) 


• 10 4 


1.99(2) 


• 10 4 


0.52 


3.82(5) 


• 10 6 


3.66(3) 


• 10 4 


3.37(3) 


• 10 4 


0.53 


1.13(2) 


•10 7 


6.36(6) 


■ 10 4 


5.87(5) 


• 10 4 


0.54 


2.97(6) 


■10 7 


1.09(1) 


■ 10 5 


1.04(2) 


• 10 5 


0.55 


9.7(3) • 


10 7 


1.95(3) 


• 10 5 


1.68(3) 


■ 10 5 


0.56 


1.9(2) • 


10 8 


3.26(6) 


■ 10 5 


3.05(9) 


• 10 5 


0.57 






5.40(8) 


■ 10 5 


4.5(1) ■ 


10 5 


0.58 






9.2(2) • 


10 5 


8.3(2) • 


10 5 


0.59 






1.61(2) 


• 10 6 


1.41(3) 


• 10 6 


0.60 






2.67(5) 


• 10 6 


2.64(8) 


• 10 6 
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One sees by comparing the results of Tables I and II 
that in most cases the two agree within 20% for temper- 
atures not too close to the critical temperature. 

A commonly used first approximation to the descrip- 
tion of the time scales of activated processes is Arrhenius' 
law, which states that the typical time scale r increases 
exponentially with the height of the (free) energy barrier 
AF = F(M = 0) - F(M = M ), with a prefactor / that 
depends in a mild way on temperature and system size, 



r = /(/?,B,L)exp(/?AF). 



(14) 



If we straightforwardly use the free energy F{M) as de- 
fined before we can obtain the heights of the free-energy 
barriers directly from the histograms of the magnetiza- 
tion distribution. To check the accuracy of this, we have 
plotted in figure 6 the ratio of r and exp((3AF) as a func- 
tion of inverse temperature /3J, for various system sizes. 
Clearly, if the prefactor / is assumed to be constant, this 
simple approximation fails to even predict the magneti- 
zation reversal times within an order of magnitude. 

A little thought reveals that, in the present case, this 
way of determining the free energy barrier is not en- 
tirely satisfactory. The reason is that for magnetization 
values around M = Mq the free energy increases with 
system size because the range of values through which 
M typically fluctuates, increases as the square root of 
system size. As a consequence the probability of find- 
ing any specific M- value decreases. Around M = on 
the other hand, no similar effect occurs, due to the wide 
plateau of the free energy as function of M. As a re- 
sult of this the decrease of probability of finding a given 
magnetization value for given relative position of the two 
interfaces is precisely compensated by the probabilities 
of finding this same magnetization value at different rel- 
ative positions. Therefore it makes more sense in this 
case to define the Arrhenius factor in an alternative way 
as exp(/3F(Af = 0)). With this definition, where again 
(5F(M) = — \n(P eq (M)), it may be interpreted indeed as 
the equilibrium probability of finding two interfaces di- 
viding the system into equal areas. But please notice that 
if, instead of a plateau, the free energy has a maximum 
of a width small compared to the range of magnetization 
fluctuations, the first definition of the Arrhenius factor is 
to be preferred. 

With the second definition the expression for the decay 
time in terms of the Arrhenius factor becomes 



/'(/3,B,L)exp(/3F(M = 0)), 



(15) 



in which again one hopes that the prefactor /' depends 
on temperature and system dimensions only in a mild 
way. On the basis of (8) and (13) we may conclude that 
for large systems /' should simply be proportional to L 
and independent of B. 

Figure 7 shows /' as a function of inverse temperature 
for several system sizes. For L = 32 and 64 both the 
approximate independence of B and the proportionality 
to L are quite well-confirmed, especially if one takes into 



account that the plateau widths for these system sizes 
are notably less than L (see figure 2). For L = 16 the 
plateau becomes so narrow that the above predictions do 
not apply. 
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FIG. 6. Ratio of the magnetization reversal times r and 
exp(/3AF), as a function of inverse temperature /3J, for vari- 
ous system sizes. The squares represent systems with 5 = 8, 
the diamonds represent systems with B = 16, and the circles 
represent systems with B — 32. 
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FIG. 7. Ratio of the magnetization reversal times r and 
exp(/3F(0)) ee P eq (M = 0) _1 , as a function of inverse tem- 
perature [3J , for various system sizes. The squares repre- 
sent systems with L = 64, the circles represent systems with 
L = 32, and the diamonds represent systems with L = 16. 
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D. most slowly decaying mode 

To obtain an estimate for the magnetization reversal 
times, we first estimated the most slowly decaying mode 
in eq. (7). Besides comparing the magnetization rever- 
sal times, we can also obtain the most slowly decaying 
cigenmode Po(M) directly from our Monte Carlo simu- 
lations by measuring the probability difference between 
the occurrence of a certain magnetization M and the oc- 
currence of the opposite magnetization — M, averaged 
over a time scale comparable to the magnetization rever- 
sal time r. Since Po decays much slower than the other 
antisymmetric (in M) modes, it will give the dominant 
contribution to this probability difference. Figure 8 com- 
pares the most slowly decaying eigenmodc as obtained 
with eq. (7) with the direct Monte Carlo measurements, 
for the 16 x 64 system at f3J = 0.5. 




500 



FIG. 8. The most slowly decaying eigenmode Po(M) di- 
vided by the equilibrium distribution exp(— f3F(M)) for the 
16 x 64 system at (3 J — 0.5, obtained in two different ways as 
described in the text. Since this function is antisymmetric in 
M, the left hand side of the figure (M < 0) has been left out 
for clarity. 

V. DISCUSSION 

Our simulations confirm the global picture we sketched 
for the process of magnetization reversals in the Ising 
model with stochastic dynamics: a large cluster of op- 
posite magnetization, originating through a fluctuation 
develops into a pair of interfaces. These interfaces dif- 
fuse around the system and annihilate, leaving the sys- 
tem in the oppositely magnetized phase. Quantitatively 
this process may be described to a good approximation 
as a diffusion process in a one-dimensional space with 
a coordinate describing the total magnetization in the 
system. 



A few remarks should be made here. 

• First of all it may seem remarkable that in almost 
all cases our theoretical prediction gives a shorter 
reversal time than the simulations. Since the the- 
ory neglects processes that enhance the reversal fre- 
quency, such as the formation of more than two 
simultaneous interfaces and spurious passages of 
M = (see below), one might expect an overesti- 
mate of the reversal time rather than an underesti- 
mate. The explanation of the latter comes primar- 
ily from that part of the process in which a growing 
cluster of opposite magnetisation transforms into a 
band around the cylinder. As the cluster will at 
first be typically circular in shape (if the tempera- 
ture is not too low) it has to deform into more ellip- 
tical shape, with a longer interface than the band, 
before the latter can be formed. In figure 2 the 
signs of this are very clear in the form of little shoul- 
ders on the sides of the plateaus in the free ener- 
gies. The transition rates Tm>,m for the formation 
of these elliptical shapes will not be proportional 
to B but much smaller, since this growth mainly 
occurs along the short sides of the cluster. The 
contributions from these terms in the denominator 
of (8) are therefore larger than estimated. Numer- 
ically they are important as long as the plateau in 
the free energy is not too wide. This then will give 
rise to a notable decrease of the reversal frequen- 
cies from our estimated values. Correcting for this 
in the theoretical expression would require better 
estimates of Tm',m in the free energy shoulders, 
which seem fairly complicated to obtain. Also the 
growth rates in the magnetization regions beyond 
the shoulders, where the growing clusters are even 
smaller, are overestimated by (13), but these are 
weighted less as the free energies there are smaller. 

• The above considerations clearly reveal the possi- 
bility that the coordinate parametrizing the "reac- 
tion path" (here the magnetization) changes non- 
monotonically along this reaction path; one could 
imagine that the growing cluster, before turning 
into a strip, typically decreases in size for a while. 
In reality this does not seem to happen, but if this 
were the case, parts of the reaction path with dif- 
ferent cluster shape but equal magnetization would 
be lumped together, and the diffusion process along 
the cluster size coordinate would no longer corre- 
spond to the actual path taken by the cluster. Ob- 
viously in such a case the reaction coordinate has 
to be redefined in such a way that the new coordi- 
nate is monotonic along the reaction path indeed. 
But in complicated situations it may not always be 
clear what is a proper choice for such a coordinate. 

• As a consequence of magnetization fluctuations in 
the bulk, passages of M = will be registered typi- 
cally already a short while before the area between 
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the interfaces reaches the value LB/ 2 correspond- 
ing to an equal division of the system between areas 
of positive and negative magnetization. The rea- 
son is that the typical time scale for fluctuations of 
the bulk magnetization is much shorter than that 
for interface diffusion, so for each location of the 
interfaces the whole range of accessible magnetiza- 
tion values typically will be scanned. Now in most 
cases this just will give rise to a negligible shift of 
the time at which the situation of equal areas (this 
is the physically relevant criterion) is reached, but 
occasionally it may happen that M = occurs, 
but the system returns to the pure state it came 
from without ever reaching equal areas. To check 
the importance of these events we may compare the 
average first passage time from — M to M to that 
from ±Mo to M = 0. If the effect of magnetization 
fluctuations is negligible the ratio of these should be 
2, otherwise it ought to be larger. Note that for the 
systems considered here the effect of magnetization 
fluctuations on the first passage time from — M to 
M is much smaller than that on the first passage 
time from —Mo to M = 0, because the probabil- 
ity of returning to —Mq without reaching the pure 
state of positive magnetization, once M = Mo has 
been reached, is extremely small. We have mea- 
sured this ratio for several system sizes and tem- 
peratures, excluding those cases where the occur- 
rence of multiple interfaces is likely, and found that 
the mean value is 1.90, with a standard deviation 
of 0.15. 

• It is of interest to investigate how the reversal times 
depend on the system size parameters. From (8) we 
may conclude that for long systems (but not so long 
that multiple interface pairs will occur frequently) 
the reversal time will become independent of the 
system length L. For this one should notice first of 
all that for M-values on the plateau exp(/3F(M)) is 
proportional to l/L, because the number of ways a 
pair of interfaces may be placed such that the aver- 
age magnetization equals M, is proportional to L. 
On the other hand the summation over M, which is 
dominated by M-values on the plateau, gives rise 
to a factor close to L. Hence to first approxima- 
tion the reversal frequency v is independent of L. 
Remarkably this independence of L in fact is ob- 
served even better by the numerical results than 
by the theoretical estimates. 

In section IV C it was noted already that the prod- 
uct /' of the reversal frequency and the Arrhcnius 
factor exp(— /3F(0)) depends on L and (3J, but is 
approximately independent of the system width B. 
Notice that square systems may be included in this 
comparison without problem, as the extra factor of 
2 due to the two possible orientations of the inter- 
faces, is properly accounted for in the Arrhenius 
factor. 



• We already indicated repeatedly that our theory 
may be applied only if the probability of having 
more than two interfaces around the cylinder may 
be neglected compared to the probability of having 
just two such interfaces. For this to be case one has 
to require Zi/Zq <C 1. For large enough systems 
this condition may be rewritten as L <C exp(/3aB), 
with a the surface tension of the interface. For large 
systems violation of this condition requires extreme 
aspect ratios, so under normal conditions it will be 
satisfied, unless the temperature is very close to the 
critical one, at which the surface tension vanishes. 

We are presently applying the methods described here 
to a study of nucleation rates in metastable states. We 
hope to report on this before long. 
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